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This  document  is  the  final  technical  report, 
prepared  under  contract  DAAK70-82-C-0011  for  the  U.S. 
Army  Engineer  Topographic  Laboratories,  Fort  Belvoir, 
VA.  The  contracting  officer's  technical  represen¬ 
tative  is  Mr.  William  A.  Allen,  ETL-TD-EA. 
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1.  INTRODUCTION 

1 . 1  BACKGROUND 

The  development  of  high  precision  self-contained  iner¬ 
tial  navigation  systems  over  the  past  two  decades  has  made 
possible  rapid  positioning  and  gravity  surveys.  Such  survey 
tools  are  used  in  both  airborne  and  land-based  configurations. 
Gravity  quantities  can  be  estimated  with  moderate  accuracy 
using  appropriate  models  and  modern  data  processing  techniques 
in  conjunction  with  the  inertial  survey  data.  Gravity  gradi- 
ometers,  operating  independently,  or  with  an  inertial  naviga¬ 
tion  system,  promise  greater  estimation  accuracy  in  comparison 
to  that  obtained  by  conventional  inertial  systems.  Even  more 
accurate,  but  slower  and  more  difficult  to  use  for  wide  area 
coverage,  is  the  astrogeodetic  methodology  for  determining 
deflections  of  the  vertical. 

Accurate  point  estimates  can  be  obtained  for  the  grav¬ 
ity  anomaly  using  gravimeters.  Because  deflections  can  be 
computed  from  fields  of  anomaly  data,  gravity  measurements  are 
also  an  approach  to  determining  deflections  of  the  vertical. 

Little  attention,  however,  has  been  given  to  the  accu¬ 
racy  potential  derived  from  processing  multisensor  data  from 
the  varied  sources  described  above.  Also,  the  accuracy  improve¬ 
ment  obtained  by  optimally  combining  multitrack  data  from  moving 
base  inertial  systems  has  yet  to  be  fully  explored.  More  inves¬ 
tigative  efforts  must  be  directed  towards  post-survey  optimal 
processing  techniques  to  extract  the  maximum  information  from 
all  measurements  made  in  a  survey  area.  These  issues  are  ad¬ 
dressed  in  this  report. 
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1.2  PURPOSE  AND  SCOPE 

An  assessment  of  the  accuracy  potential  using  high 
precision  inertial  systems  with  or  without  gravity  gradiometers 
to  estimate  gravity  disturbance  quantities  was  performed  in 
Ref.  1.  In  Ref.  2,  the  use  of  astrogeodetic  or  gravimetric 
data,  in  a  single  sensor,  multimeasurement  context,  to  estimate 
deflections  of  the  vertical,  was  considered.  Documenting  the 
second  phase  of  that  effort,  this  report  explores  the  accuracy 
realizable  by  optimally  processing  astrogeodetic  data  alon^  or 
in  combination  with  gravimetric  data.  In  addition,  the  qt-  ity 
of  deflection  estimates  which  are  achievable  by  combining  Iti- 
track  data  from  inertial  and  gradiometric  systems  with  ast 
geodetic  data,  is  analyzed.  Tradeoffs  between  accuracy  a 
able  and  survey  effort  (e.g. ,  number  of  measurements)  are  made. 
Decision  making  parameters  are  identified,  with  emphasis  on 
accuracy,  cost  and  development  time.  Summaries  of  these  findings 
are  presented  in  Tables  4-1  and  4-2.  The  quantities  of  interest 
are  deflections  of  the  vertical  and  rms  deflection  estimation 
errors. 


Deflection  of  the  vertical  recovery  accuracies  are 
reported  for  two  types  of  survey  scenarios, 

•  Astrogravimetric  Survey 

•  Astrogeodetic/Moving  Base  Inertial  Sys¬ 
tems  Survey. 

Astrogeodetic  measurements  combined  with  gravimetric  measure¬ 
ments  in  a  defined  area  make  up  the  astrogravimetric  survey. 

In  the  second  survey  type,  multitrack  data  from  an  inertial 
system  and  astrogeodetic  measurements  are  combined  to  provide 
improved  deflection  estimates  between  astro  stations.  Although 
the  multitrack  aspect  of  the  inertial  survey  investigations  was 
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beyond  the  original  scope  which  TASC  had  intended  to  pursue  in 
this  study,  it  became  evident  that  meaningful  comparisons  amonj 
survey  alternatives  required  that  inertial  surveys  be  treated 
in  the  same  context  as  those  providing  astrogeodetic  and  gravi¬ 
metric  data.  This  context  involves  a  rigorous  treatment  of 
minimum  variance  processing  techniques  applied  to  data  taken 
in  two-dimensions. 


Three  types  of  moving  base  inertial  systems  are  stud¬ 
ied  in  this  effort.  They  are:  a  conventional  inertial  system 
(the  baseline  of  Ref.  1,  representative  of  the  Rapid  Gravity 
Survey  System  [RGSS]),  a  star  tracking  inertial  system  (the 
Northrop  NAS-26  astroinertial  navigator) ,  and  a  moving-base 
gravity  gradiometer  system.  The  use  of  the  Ref.  1  baseline 
scenario,  a  helicopter  traverse  over  a  straight,  35  km  track, 
is  maintained  in  this  portion  of  the  study  for  those  cases 
involving  survey  tracks.  Thus  comparisons  with  the  previous 
single  track  results  are  most  meaningful.  Note  that  actual 
data  is  not  processed  or  analyzed  in  this  report.  Only  the 
estimation  accuracy  potential  of  the  methods  for  the  defined 
survey  arrangements  is  determined. 


1.3  OVERVIEW  OF  REPORT 

This  report  is  organized  as  follows:  Chapter  2  dis¬ 
cusses  the  deflection  estimation  accuracy  associated  with  astro- 
gravimetric  surveys.  Chapter  3  addresser  the  astrogeodetic/ 
moving  base  survey  configurations.  Costs  and  logistical  issues 
are  considered  in  Chapter  4  which  also  extends  to  a  summary  of 
important  single  sensor  technological  alternatives,  for  deter¬ 
mining  deflections  of  the  vertical.  A  summary  is  provided  in 
Chapter  5. 
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2.  ASTROGRAVIMETRIC  SURVEYS 

2.1  INTRODUCTION 

In  this  chapter  the  residual  error  in  the  estimate  of 
the  deflection  of  the  vertical  at  a  point,  given  measurements 
of  both  the  deflection  and  gravity  anomaly  in  the  vicinity  of 
the  point  and  an  appropriate  worldwide  gravity  data  base,  is 
determined.  This  is  performed  by  simulating  the  errors  in 
measuring  the  various  astrogeodetic  and  gravimetric  quantities. 
In  the  actual  processing  algorithm,  the  measurements,  consist¬ 
ing  of  astrogeodetic  and  gravimetric  data  taken  over  a  speci¬ 
fied  survey  area,  are  optimally  combined  to  yield  the  desired 
deflection  estimates.  Since  anomaly  and  deflection  are  corre¬ 
lated  to  some  degree  over  the  survey  area,  the  combination  of 
both  types  of  measurements  results  in  final  error  residuals 
which  are  smaller  than  those  derived  from  either  measurement 
alone.  The  numerical  results  presented  here  quantify  the  degree 
of  improvement  obtained  by  mixing  astrogeodetic  and  gravimetric 
measurements  in  an  optimal  scheme  and  define  the  survey  effort 
required  to  achieve  a  specific  accuracy  goal.  The  sensitivity 
of  the  residual  errors  to  such  factors  as  survey  area  geometry 
and  gravity  field  characteristics  is  determined  and  provides 
insight  into  the  advantages  of  this  survey  strategy. 

In  the  next  section  the  simulation  models  used  in  the 
analysis  are  presented  and  described.  These  include  the  grav¬ 
ity  field  models  and  the  sensor  (measurement)  error  models.  A 
description  of  the  survey  area  geometry  is  given  in  Section  2.3 
and  the  methodology  employed  in  determining  the  error  residuals 
is  detailed  in  Section  2.4.  Numerical  results  of  this  portion 
of  the  study  are  presented  in  Section  2.5. 
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2.2  SIMULATION  MODELS 

The  mathematical  models  developed  in  this  report  are 
statistical  descriptions  of  random  processes.  As  such  they 
are  used  to  describe  the  ensemble  or  average  behavior  of  the 
quantities  modeled.  Auto-  and  cross-correlation  matrices  used 
in  the  numerical  formulation  are  computed  from  analytic  ex¬ 
pressions  for  the  correlation  functions.  Frequency  domain 
descriptions  (power  spectral  densities)  are  also  employed. 

2.2.1  Gravity  Disturbance  Models 

The  third-order  Markov  analytic  gravity  model  described 
and  used  in  Ref.  1  (see  also  Ref.  4)  is  retained  for  the  Phase  II 
studies  to  characterize  the  gravity  disturbance  field  (raw 
field).  However,  for  the  astrogravimetric  portion  of  the  studies, 
the  model  parameters  have  been  modified  to  account  for  distant 
zone  data  as  described  below. 

It  is  assumed  that  appropriate  data  from  the  DMA  gravity 
library  data  base  is  to  be  used  for  "distant  zone"  coverage 
when  processing  gravimetric  surveys.  This  data  is  available 
as  area  means  in  one  deg  x  one  deg  blocks  all  over  the  world, 
and  at  smaller  intervals  in  some  specific  areas.  In  ocean 
areas,  the  DMAAC  data  base  includes  the  satellite  altimetery 
information  provided  by  GEOS- 3  and  SEASAT-1.  For  the  (present) 
purpose  of  performing  a  survey  effort  and  estimation  accuracy 
assessment,  the  statistics  of  the  anomaly  errors  are  modeled. 

These  are  based  on  tabulated  uncertainty  estimates  which  are 
maintained  as  part  of  the  gravity  data  base.  In  Ref.  3,  numer¬ 
ical  models  are  developed  to  represent  the  power  spectral  den¬ 
sity  of  the  free-air  anomaly  uncertainty.  Separate  representa- 


The  model  is  "Markov",  i.e.,  single  step  predictive/recursive, 
for  processes  modeled  along  a  single  track. 
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tions  are  developed  for  (1)  the  entire  data  base,  yielding  a 
world-wide  (WW)  model,  and  (2)  for  data  in  and  near  the 
Continental  U.S.  (CONUS)  only,  a  CONUS  model. 

The  self-consistent,  third-order  Markov  formulation 
of  the  WW  and  CONUS  models  also  provides  some  important  compu¬ 
tational  advantages  over  the  numerical  models  of  Ref.  3,  the 
major  one  being  that  the  auto-  and  cross-correlation  for  deflec¬ 
tions  of  the  vertical  and  anomaly  are  analytically  computed 
(as  a  function  of  geometry)  from  simple  formulas  tabulated  in 
the  cited  references.  For  the  astrogravimetric  survey  the  acf 
(autocorrelation  function)  form  of  the  model  is  applicable,  as 
opposed  to  the  state-space  form  needed  for  temporal  models 
(e.g.,  inertial  navigation  system  simulation). 


The  anomaly  power  spectral  density  (PSD)  function  for 
the  third-order  Markov  model  is  given  by 


<u0 


15n  p3  o2  u)2 

7 zr  -  277 72 

(P  +  U)  ) 


(2.2-1) 


where , 

uj  ~  two-dimensional  radial  spatial  frequency 
( rad/km) 

p  ~  corner  frequency  (inverse  of  character¬ 
istic  distance) 

a  ~  rms  value  of  the  gravity  anomaly 

o 

The  two-dimensional  PSD  given  by  Eq.  2.2-1  is  needed  since 
gravimetric  uncertainty  in  the  entire  horizontal  plane  is  to 
be  determined.  The  "corner  frequency",  or  its  inverse,  the 
characteristic  distance,  divides  the  gravity  anomaly  spectrum 

•fa 

Self-consistency  is  discussed  at  length  in  Ref.  U. 
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into  two  regions.  Most  of  the  anomaly  energy  occurs  at  fre¬ 
quencies,  u><g.  Anomaly  frequencies  greater  than  g  contain 
decreasing  amounts  of  energy.  Hence,  a  long  characteristic 
distance  (low  corner  frequency)  implies  a  gravity  disturbance 
field  which  varies  slightly  over  short  distances  while  a  short 
characteristic  distance  (high  corner  frequency)  implies  a  grav¬ 
ity  field  which  changes  rapidly  over  short  distances.  The 
former  characteristics  are  more  desirable  from  an  estimation 
viewpoint,  since,  for  a  gravity  disturbance  field  which  varies 
only  slightly  over  short  distances,  an  estimate  of  a  gravity 
quantity  at  a  point,  given  nearby  measurements,  can  be  made 
with  reasonable  certainty.  However,  for  a  field  which  varies 
significantly  over  short  distances,  this  uncertainty  is  larger, 
since  less  information  can  be  obtained  from  the  measurements. 

The  post-measurement  estimation  uncertainty  approaches  that  of 
the  reference  field  for  very  short  anomaly  characteristic 
distances . 

The  parameters  of  the  third-order  Markov  model  can  be 
selected  to  represent  either  the  worldwide  or  the  CONUS  grav¬ 
ity  anomaly  models  of  Ref.  3.  These  parameter  values  are  listed 
in  Table  2.2-1. 

TABLE  2.2-1 

GRAVITY  DATA  BASE  UNCERTAINTY  MODEL  PARAMETERS 


PARAMETER 

VALUE 

WW  MODEL 

CONUS  MODEL 

og  (mgal ) 

38 

14 

g"1  (km) 

35 

16 

*An  extreme  version  of  this  situation  occurs  in  mountainous 
terrain.  Even  when  dense  gravity  data  is  available,  adequate 
treatment  requires  additional  information  in  the  form  of 
topography  models.  Modern  thinking  in  this  regard  is  set 
forth  in  Ref.  7. 
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2.2.2  Sensor  Error  Models 


The  data  processing  technique,  described  in  Section  2.3, 
requires  statistical  descriptions  (covariance  matrices)  of  the 
measurement  errors  to  compute  the  estimates  and  their  error 
statistics.  Error  models  are  defined  for  the  astrogeodetic 
and  gravimetric  measurements  which  comprise  the  astrogravimet- 
ric  survey.  Simple  models  have  been  chosen,  but  the  parameter 
values  are  representative  of  actual  measurement  processes 
(Ref.  2).  These  are  theodolite  measurements  and  gravimeter 
measurements  reduced  to  the  geoid. 

A  first-order  Markov  noise  model  is  used  to  describe 
the  astrogeodetic  measurement  errors.  The  autocorrelation 
function  of  this  process  is  given  by 

4>(x)  =  o2  e"x/d  (2.2-2) 

2 

where  a  is  the  variance,  x  is  the  shift  distance,  and  d  is 
the  correlation  distance.  This  model  accounts  for  spatial 
correlations  between  astro  measurements  at  different  sites. 

Gravity  anomaly  measurement  errors  are  modeled  using 
a  white  noise  description.  The  appropriate  random  process  is 
described  by 


<t> (x)  =  A  6 (x)  (2.2-3) 

where  A  is  the  strength  of  the  process,  x  is  the  shift  distance, 
and  6  is  the  Dirac  impulse  function. 


*An  example  of  one  correlation  mechanism  is  anomalous  atmospheric 
refraction.  Weather  cells  or  local  conditions  can  cause  simi¬ 
lar  refraction- induced  errors  over  areas  of  one  km  or  more. 
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A  summary  of  the  instrument  (measurement)  error  models 
is  provided  in  Table  2.2-2.  Values  for  the  error  model  param¬ 
eters  used  in  the  present  study  are  also  listed.  The  anomaly 
error  includes  the  effects  of  reduction  and  downward  continua¬ 
tion  of  gravimeter  measurements  to  the  geoid. 


TABLE  2.2-2 
SENSOR  ERROR  MODELS 


SENSOR 

(MEASUREMENT) 

MODEL 

PARAMETER  VALUES 

Astrogeodetic 

First-Order 

0 . 3  sec  rms 

Deflection  of 
the  Vertical 

Markov 

correlation  distance  =  1  km 

Gravimeter 

Uncorrelated 

1  mgal 

(Anomaly) 

post-survey 

reduction 

errors 

2.3  SURVEY  GEOMETRY 

! 

The  astrogravimetric  survey  geometry,  shown  in 
1  Fig.  2.3-1,  is  a  square  pattern  of  measurement  points  in  the 

horizontal  plane.  It  represents  either  a  very  local  survey  or 
|  a  "cell"  which  can  be  repeated  over  a  large  area.  Astrogeo- 

!  detic  measurements  (North  and  East)  are  made  at  the  corners  of 

the  cell  and  may  or  may  not  be  made  at  the  cell  midpoint. 
Deflection  of  the  vertical  survey  error  is  evaluated  at  the 
center  of  the  cell.  The  simulation  studies  address  an  arbi¬ 
trary  number  of  gravity  anomaly  measurements  distributed  over 
the  area  of  the  cell;  nine  are  illustrated  in  Fig.  2.3-1.  The 
number  of  anomaly  measurements  (e.g.,  the  densif ication  of 
anomaly  data)  is  varied  in  the  analysis  to  observe  the  change 
in  the  midpoint  deflection  of  the  vertical  residual  error  for 
a  given  astro-measurement  spacing. 
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Figure  2.3-1  Typical  Survey  Geometry 


2.4  ACCURACY  EVALUATION  METHODOLOGY 

The  analysis  herein  is  based  on  the  use  of  optimal 
(minimum  variance)  data  processing  techniques  to  estimate  grav¬ 
ity  quantities  from  zero  mean  field  measurements.  In  this 
respect  a  vector  quantity  x  is  estimated  from  a  vector  of 
measurements  z  by 


x 


=  4> 


xz 


(2.4-1) 


where 


=  cross-covariance  between  x  and  z 

AZ 

*  =  covariance  of  z 

zz  — 

The  covariance  of  the  associated  residual  estimation  errors  is 
then  computed  from 
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(b  •  rt)  -  A  (b  A 
Yxx  Yxx  Yxz  YZZYZX 

where  is  the  a  priori  covariance  of  x. 

If  measurements  (z)  of  x  are  made,  then 


(2.4-2) 


d)  =  d»  +  R 

YZZ  YXX  X 


and  4>  =  <ft 


xz 


xx 


where  Rx  is  the  measurement  error  covariance  matrix.  This 
processing  technique,  also  known  as  spatial  collocation  (see 
Ref.  5),  forms  the  basis  for  the  accuracy  analysis  performed 
here.  Its  main  advantages  are: 


•  Minimum  variance  estimates  are  determined 

•  No  restrictions  are  placed  on  the  spacing 
or  locations  of  measurement  or  estimation 
points 

•  Systematically  incorporates  multisensor 
data 

•  Computationally  efficient  for  a  moderate 
number  of  measurements. 


A  more  elaborate  extension  of  this  formulation  involves  com¬ 
puting  the  Fourier  transforms  of  the  quantities  in  Eq.  2.4-2 
and  performing  the  computations  in  the  frequency  domain.  In 
this  way  the  measurement  covariance  matrix  to  be  inverted  becomes 
diagonal,  thereby  streamlining  the  computations  for  large  numbers 
of  measurements.  Such  a  method  is  described  in  Ref.  6.  For 
the  purposes  of  this  study,  however,  spatial  collocation  is 
adequate . 


The  covariance  matrices  of  Eq.  2.4-2  are  analytically 
computed  from  the  third-order  Markov  expressions  tabulated  in 
Ref.  4  for  deflections  of  the  vertical  and  gravity  anomaly. 
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As  previously  discussed,  the  free  parameters  are  chosen  to  fit 
the  Markov  model  to  the  specified  reference  models  of  Ref.  2 
(world-wide  or  CONUS).  The  survey  geometry  defines  the  vector 
distances  between  measurement  and  estimation  points,  and  these 
distances  are  used  to  construct  the  covariance  matrices.  The 
measurement  error  covariance  matrices  are  computed  from  the 
error  models  of  Section  2.2.2. 

Equation  2.4-2  may  be  solved  once  for  all  of  the  meas¬ 
urements,  astrogeodetic  and  gravimetric,  in  a  batch  formula¬ 
tion.  However,  it  is  more  informative  to  use  the  sequential 
collocation  formulation  of  Ref.  5  to  process  the  data  in  a 
sensor  by  sensor  manner.  In  this  way  the  accuracy  effects  of 
one  measurement  type,  say,  of  anomaly,  can  be  seen  with  respect 
to  the  impact  of  other  measurement  types  (i.e.,  East  and  North 
deflections  of  the  vertical).  The  sequential  collocation  for¬ 
mulas  for  three  measurements  are  given  in  Appendix  A. 

I 

2.5  NUMERICAL  RESULTS 

!  I 

The  accuracy  evaluation  methodology  outlined  in  Sec- 
|  tion  2.4  is  exercised  here  to  quantify  the  performance  potential 

of  the  application  of  the  astrogravimetric  method  using  minimum 
|  |  variance  estimation.  The  simulation  models  of  Section  2.2  are 

I  employed  for  a  number  of  survey  geometry  variations  and  measure- 

■  -  ment  schemes  to  determine  the  survey  conditions  required  to 

[  achieve  a  specific  accuracy  goal.  In  this  study,  the  survey 

cell  midpoint  serves  as  the  accuracy  evaluation  point,  and,  for 
|  .  making  comparisons  and  tradeoffs,  a  baseline  accuracy  goal  of 

0.2  sec  for  the  deflection  estimate  is  established. 

It  is  remarked  in  Section  2.2.1  that  the  estimation 
accuracy  of  a  geopotential  quantity,  given  nearby  measurements, 
is  particularly  sensitive  to  the  high  frequency  content  of  the 
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gravity  disturbance  field.  This  is  shown  in  Fig.  2.5-1  where 

the  deflection  estimation  uncertainty  at  the  cell  midpoint  is 

plotted  as  a  function  of  normalized  characteristic  distance. 

The  normalization  used  is  the  characteristic  distance  of  the 

worldwide  gravity  disturbance  data  base  model.  Note  that,  in 

Fig.  2.5-1,  the  world  wide  gravity  disturbance  characteristic 

distance,  8~^,  is  taken  from  Table  2.2-1. 
ww 

Increasing  the  model  characteristic  distance  is  math¬ 
ematically  equivalent  to  moving  the  measurements  closer  together. 
For  example,  if  t  were  decreased  to  5  km,  the  appropriate  ab¬ 
scissa  in  Fig.  2.5-1  would  be  2.  In  the  limit  of  zero  spacing 
(or  equivalently,  infinite  characteristic  distance)  the  error 
quantity  is  a  bias  and  the  residual  error  of  the  estimate  is 
simply  the  (uncorrelated)  measurement  error  divided  by  the 
square  root  of  the  number  of  measurements.  In  the  range  of 
cases  described  by  Fig.  2.5-1,  the  measurement  error  for  each 
of  the  deflection  measurements  is  0.3  sec  and  the  astro- 
measurement  error  correlation  distance  is  small  compared  to 
the  measurement  spacing  (the  astro-measurement  errors  are  ef¬ 
fectively  uncorrelated  with  each  other).  Therefore  all  of  the 
curves  of  Fig.  2.5-1  approach  0.15  sec  as  the  gravity  model 
characteristic  distance  increases. 

The  number  of  anomaly  measurements  involved  is  also  a 
simulation  parameter.  As  indicated  in  Fig.  2.5-1,  the  number 
of  anomaly  measurements  is  important  in  reducing  the  deflection 
estimation  uncertainty  when  the  gravity  disturbance  field  con¬ 
tains  significant  short  wavelength  energy.  For  gravity  distur¬ 
bance  fields  with  modest  short  wavelength  content,  nearby  anomaly 
measurements  are  of  little  benefit.  A  deflection  estimate  re¬ 
ceives  maximum  information  from  an  anomaly  measurement  when  the 
two  points  are  separated  by  about  one  characteristic  distance. 
This  aspect  is  typified  by  the  small  local  minimum  in  the  bottom 
curve  of  Fig.  2.5-1. 
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A  N.  E  DEFLECTION  MEASUREMENT 
X  ESTIMATE  POINT 


ANOMALY  MEASUREMENTS 
IN  AREA  (STUDY  VARIABLE) 


GRAVITY  DATA  BASE  MODEL  CHARACTERISTICS  DISTANCE 


Figure  2.5-1  Sensitivity  of  Deflection  of  the  Vertical 

Estimate  Uncertainty  to  Gravity  Disturbance 
Field  Characteristic  Distance 


An  important  factor  in  determining  astrogravimetric 
survey  effort  requirements  involves  the  tradeoff  between  the 
deflection  measurement  spacing  and  gravimeter  measurement  den- 
sification  within  the  survey  area.  The  variations  of  the  mid¬ 
point  North  deflection  residual  error  (rms)  with  these  factors 
are  summarized  in  Fig.  2.5-2  using  the  worldwide  gravity  data 
base  uncertainty  model  with  no  astrogeodetic  measurements  at 
the  estimate  point.  It  is  clearly  seen  that  a  larger  astro 
data  spacing  can  be  tolerated  at  the  expense  of  increased  den- 
sification  of  the  gravimetric  measurements.  No  anomaly  data 
is  required  to  achieve  the  estimation  accuracy  goal  if  the 
astro-measurement  spacing  is  small.  As  the  astro  measurement 
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Figure  2.5-2  Alternative  Astrogravimetric  Measurement 

Strategies  without  Midpoint  Deflection 
Measurement 

separation  distance  increases,  an  increasing  number  of  anomaly 
measurements  are  required  to  achieve  the  goal.  If  the  astro- 
measurement  spacing  is  large,  then  no  reasonable  number  of 
anomaly  measurements  will  result  in  the  specified  accuracy 
level.  Fewer  gravimetric  measurements  are  required,  or  alter¬ 
natively,  larger  astro-measurement  spacings  can  be  used  if 
deflection  measurements  are  made  directly  at  the  estimation 
point.  This  is  shown  in  Fig.  2.5-3. 

Similar  graphs  are  presented  in  Figs.  2.5-4  and  2.5-5 
for  the  CONUS  gravity  reference  field.  The  results  illustrated 
in  Figs.  2.5-4  and  2.5-5  are  indicative  of  typical  astrogravi¬ 
metric  survey  performance  to  be  expected  on  the  U.S.  land  mass. 
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e  2.5-3  Alternative  Astrogravimetric  Measurement 
Strategies  with  Midpoint  Deflection 
Measurement 


A  K.  C  OEFLECTKMt  MEASUREMENTS 
X  ESTIMATE  POINT 


ANOMALY  MEASUREMENTS 
M  AREA  (STUDY  VARIABLE) 


Figure  2.5-4  Alternative  Astrogravimetric  Measurement 

Strategies  without  Midpoint  Deflection 
Measurement 
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Figure  2.5-5  Alternative  Astrogravimetric  Measurement 

Strategies  with  Midpoint  Deflection 
Measurement 


For  a  given  set  of  survey  parameters,  the  CONUS  area  and  data 
base  results  in  estimation  accuracies  similar  to  those  derived 
from  the  worldwide  field.  This  is  not  surprising  since  the 
measurement  densif ication  of  the  survey  models  studied  here 
greatly  exceeds  the  average  densif ication  of  existing  world¬ 
wide  or  CONUS  data. 
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ASTR0GE0DETIC/M0V1NG  BASE  INERTIAL  SYSTEM  SURVEYS 

3 . 1  INTRODUCTION 

In  this  chapter  the  deflection  estimation  accuracy 
achievable  by  optimally  combining  data  from  inertial  systems 
and  astrogeodetic  measurements  is  explored.  A  major  output  of 
this  effort  is  a  formulation  to  combine  mul ti track  data  from 
moving  base  inertial  systems,  along  with  the  astrogeodetic 
measurements,  to  "re-estimate"  deflections  of  the  vertical  on 
the  separate  survey  tracks  in  a  post  survey  regional  adjustment. 
The  astrogeodetic  measurements  serve  as  endpoint  calibration 
values  for  the  baseline  helicopter  traverse  studied  herein. 

This  work  is  a  direct  continuation  of  the  covariance  studies 
made  in  Ref.  1  for  a  single  traverse  of  an  inertial  surveying 
configuration.  The  baseline  error  models  (gravity,  sensor, 
etc.)  and  the  survey  parameters  of  Ref.  1  are  retained  for  use 
here . 

Three  types  of  moving  base  inertial  systems  are  simu¬ 
lated  for  making  accuracy  potential  comparisons.  They  are 

1)  a  conventional  inertial  system  typified 
by  the  RGSS  (the  baseline  of  Ref.  1) 

2)  a  gradiometer-aided  inertial  system  (the 
Bell  configuration  is  modeled) 

3)  the  Northrop  NAS-26  astroinertial 
navigation  system. 

The  NAS-26  utilizes  a  star- tracker ,  integral  with  the  navigator, 
for  making  frequent  star  sightings  to  calibrate  the  gyros  and 
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bound  the  position  error.  Although  not  configured  for  gravity 
surveys,  the  NAS-26  hardware  is  currently  available.  Some 
software  and  a  zero  velocity  update  capability  would  have  to 
be  added  for  survey  use.  The  gradiometer-aided  system  is  not 
presently  available  for  use  in  gravity  surveys,  but  the  abil¬ 
ity  of  a  gradiometer  configuration  to  estimate  deflections  of 
the  vertical  is  assessed  for  future  survey  planning. 


3.2  SIMULATION  MODELS 

The  gravity  disturbance  field  model,  the  navigation 
error  model,  and  the  sensor  error  models  used  here  are  iden¬ 
tical  to  those  of  Ref.  1  wherein  single  track  results  are  re¬ 
ported.  Major  characteristics  of  these  models  are  listed  here 
for  completeness  and  reference.  Also,  the  error  model  for  the 
NAS-26  is  described  in  the  following. 

3.2.1  Gravity  Disturbance  Model 

The  self-consistent  third-order  Markov  model  described 
in  the  previous  chapter  is  utilized  here  to  describe  the  statis¬ 
tics  of  the  gravity  disturbance  field.  However,  as  stated  in 
Ref.  1,  the  parameters  of  the  model  have  been  adjusted  to  fit 
the  White  Sands  gravity  anomaly  data  to  provide  a  gravity  model 
that  is  representative  of  that  region.  The  parameters,  rms 
gravity  anomaly  and  characteristic  distance,  are  presented  in 
Table  3.2-1.  The  Markov  gravity  disturbance  model  is  used  in 
its  velocity-invariant  (any  vehicle  velocity  can  be  simulated) 
state-space  form  for  covariance  analysis,  and  contains  the 
structure  necessary  to  formulate  the  gravity  gradients  needed 
to  model  the  gradiometer  measurements.  Details  of  the  gravity 
model  can  be  found  in  Appendix  B  of  Ref.  1. 
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TABLE  3.2-1 

WHITE  SANDS  GRAVITY  MODEL  PARAMETERS 


PARAMETER 

VALUE 

UNITS 

Anomaly  rms ,  a 

Characteristic  distance,  p  ^ 

15.8 

17.9 

mgal 

km 

3.2.2  Navigation  System  Error  Model 

The  main  components  of  an  inertial  surveying  system 
are  gyros,  accelerometers,  a  stabilized  platform,  support  elec 
tronics,  power  supplies,  environmental  protection,  and  the 
appropriate  software  to  compute,  using  sensed  accelerations, 
the  position  and  velocity  of  the  carrying  vehicle.  The  error 
equations  for  a  full  three-axis  unaided  inertial  navigation 
system  are  given  below.  The  position,  velocity,  and  misalign¬ 
ment  errors,  coordinatized  in  the  local-level  frame  (North, 
East,  Down)  are  driven  by  accelerometer  errors,  gravity  dis- 


turbances ,  and  gyro  drifts. 

fir  =  fiv  -  u)ri  x  fir 
—  —  —EL,  — 

(3.2-1) 

fiv  =  u-fig~4)xA  -  (wjL  +  ft)  x  fiv 

(3.2-2) 

tji  =  -  x  +  € 

(3.2-3) 

where 


fir  =  position  error  (computed-minus- true ) 
fiv  =  velocity  error  (computed-minus-true) 
ij>  =  misalignment  (platform- to-computer) 
fig  =  gravity  error  (true-minus-reference) 
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A  =  specific  force;  i.e.,  vehicle  acceleration  due  to 
all  forces  acting  except  gravity  (ideal  acceler¬ 
ometer  output) 

Q  =  angular  rate  of  earth- fixed  axes  with  respect  to 
inertial  space 

u  =  accelerometer  error 

£  =  gyro  drift 

-IL  =  an£ular  rate  of  the  local-level  frame  (L)  with 
respect  to  inertial  space  (I) 

Wp,  =  angular  rate  of  the  local-level  frame  with  respect 
to  earth  (E) 

Endpoint  calibration  accuracies  for  the  survey  traverse 
simulations  are  listed  in  Table  3.2-2.  These  accuracies,  used 
at  both  the  beginning  and  the  end  of  the  traverse  simulation, 
model  optimal  endpoint  calibration  procedures.  Position  errors 
are  stated  relative  to  a  local  benchmark.  Further  discussion 
of  endpoint  initialization  and  termination  values  is  given  in 
Ref.  1. 


TABLE  3.2-2 

ENDPOINT  CALIBRATION  DATA 


PARAMETER 

VALUE 

UNITS 

Deflection  of  the  Vertical 

0.2 

sec 

Endpoint  Gravity  Acceleration 

0.1 

mgal 

Relative  Position  Error 

1.0 

cm 

Zero  Velocity  Error 

0.002 

km/hr 
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The  state-space  form  of  the  error  equations  allows 
straightforward  measurement  and  update  implementations  for 
zero  velocity  updates  (ZUPTS),  gradioraeter  measurements,  the 
star  sightings,  and  endpoint  calibration  procedures.  The  ex¬ 
panded  (scalar)  form  of  the  navigation  error  equations  are 
presented  in  Appendix  A  of  Ref.  1. 

3.2.3  Sensor  Error  Models 

Gyroscopes  and  accelerometers  are  the  major  components 
of  an  inertial  positioning  system  and  contribute  random  errors 
to  the  quantities  estimated.  An  error  model  for  the  inertial 
sensors  is  used  in  this  study  and  integrated  with  the  system 
error  model. 

Baseline  Sensors  -  An  initial  random  bias  drift  and  a 
white  noise  drift  are  used  to  model  the  gyro  errors.  The  param 
eters  selected  for  these  models  reflect  the  quality  of  present- 
day  production  gyros.  Values  are  given  in  Table  3.2-3.  Discus 
sion  is  provided  in  Ref.  1. 

First-order  Markov  models  are  employed  to  simulate 
accelerometer  bias  instability  and  scale  factor  errors.  A 
white  noise  accelerometer  error  is  also  included.  The  Markov 
terms  account  for  the  moderately  rapid,  yet  bounded,  error 
growth.  White  noise  adequately  models  platform  jitter,  loop 
rebalance  noise,  and  other  high-frequency  errors. 

For  the  survey  system  employing  a  gravity  gradiometer 
triad,  a  first-order  Markov  and  a  white  noise  error  model  are 
used  to  represent  the  self-noise  effects  of  each  of  the  gradi- 
ometers.  Parameter  values  of  the  models  may  be  chosen  to  simu¬ 
late  either  the  Bell  or  Draper  gradiometer  configurations. 
Models  for  the  Bell  system  are  documented  here.  They  typify 
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TABLE  3.2-3 

BASELINE  SENSOR  ERROR  MODELS 


SENSOR 

MODEL  DESCRIPTION 

VALUE 

UNITS 

Gyroscope 

Bias  Drift 

0.003 

deg/hr 

White  Noise  Drift 

0.003 

deg/hr 

(after  one 

hour) 

Bias  Instability-Markov 

Accelerometer 

Initial  Value 

2 

M8 

Time  Constant 

0.5 

hr 

Scale  Factor-Markov 

Initial  Value 

1.5 

ppm 

Time  Constant 

0.5 

hr 

White  Noise 

1.0 

Mg 

Gradiometer 

Self-noise-Markov 

(Bell) 

Initial  Value 

49 

Eotvos  (E) 

Time  Constant 

2105 

hr 

White  Noise 

86 

E2/Hz 

the  survey  results  to  be  expected  with  gradiometers  (little 
difference  is  observed  in  deflection  of  the  vertical  estima¬ 
tion  accuracy  between  the  Bell  or  Draper  systems  for  the  sce¬ 
narios  described  in  Ref.  1).  A  summary  of  the  baseline  sensor 
error  models  is  presented  in  Table  3.2-3. 

NAS-26  Astroinertial  Navigation  System  --  The  NAS -26 
astroinertial  navigation  system  consists  of  an  Inertial  Measure¬ 
ment  Unit  (IMU)  that  includes  an  enclosed  three-gimbal  stable 
reference  platform  and  a  star  tracker.  The  IMU  is  made  up  of 
two  Kearfott  Gyroflex  Mod  2  two-degree-of- freedom  dry- tuned 
rotor  gyros  and  three  Kearfott  Model  2401  single-axis  acceler¬ 
ometers.  The  entire  platform  assembly  can  revolve  about  a 
vertical  axis  in  the  base  of  the  pitch  gimbal.  The  star  tracker 
is  mounted  on  the  three-gimbal  inertial  platform  in  such  a 
manner  as  to  allow  rotation  both  in  azimuth  about  a  vertical 
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axis  relative  to  the  platform,  and  also  in  elevation  about  a 
horizontal  axis.  The  instantaneous  f ield-of-view  of  the  tele¬ 
scope  is  AO  sec . 

Astroinertial  operation  involves  using  frequent  star 
sightings  as  navigation  references.  The  ephemerides  for  61 
stars  whose  magnitudes  range  from  -1.46  to  +3.50  are  stored  in 
computer  memory  (and  updated  annually).  Assuming  cloud  cover 
does  not  pose  a  problem,  an  average  of  three  stars  per  minute 
may  be  tracked.  Normally,  six  to  ten  stars  are  available  at 
any  given  time  and  a  minimum  of  two  stars  are  available  anywhere/ 
anytime.  An  onboard  Kalman  filter  processes  the  stellar  angular 
measurements  to  estimate  updates  for  the  states  included  in 
the  onboard  filter  model  (e.g.,  position,  velocity,  tilt,  head¬ 
ing,  gyro  drift,  accelerometer  bias,  star  tracker  elevation 
bias).  Thus,  these  stellar  updates  tend  to  bound  navigation 
position  errors  and  provide  gyro  drift  rate  compensation, 
resulting  in  signficantly  improved  inertial  performance.  For 
the  study  involving  the  NAS-26  navigation  system,  error  models 
appropriate  for  the  component  instruments  are  used.  The  error 
model  parameters  for  these  sensors  are  shown  in  Table  3.2-4. 

Note  that  these  models  indicate  sensor  components  of  lower 
quality  than  the  baseline  inertial  system  discussed  in  the 
previous  section.  The  star  tracker  is  modeled  at  a  sighting 
accuracy  of  one  sec  rms  per  axis  per  measurement. 


3.3  ACCURACY  EVALUATION  METHODOLOGY 

In  this  section  the  estimation  and  the  estimation 
error  covariance  equations  are  formulated  for  multitrack  gravity 
data  processing.  The  method,  which  is  based  on  Gauss-Markov 
theory  and  spatial  collocation,  handles  the  multisensor/multi¬ 
track  data  (to  be  defined  shortly)  arrangements  systematically, 
and  ensures  minimum  variance  estimates  with  respect  to  the 
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TABLE  3.2-4 

NAS-26  SENSOR  ERROR  MODELS 


SENSOR 

MODEL  DESCRIPTION 

VALUE 

UNITS 

Gyroscope 

Initial  Drift  Rate 

0.05 

deg/hr 

Drift  Rate  Noise  (after  1  hr) 

0.05 

deg/hr 

Accelerometer 

Bias  Instability-Markov 

Initial  Value 

5 

Mg 

Time  Constant 

24 

hr 

Scale  Factor  (Bias) 

3 

ppm 

White  Noise  (10  sec  average) 

1 

Mg 

Star  Tracker 

White  Noise 

1 

sec 

data  and  the  error  models  used.  The  error  equations  determine 
the  accuracy  improvement  obtained  in  the  single-track  deflec¬ 
tion  estimates  by  employing  the  new  formulation  in  a  post- survey 
regional  adjustment.  The  formulation  also  accounts  for  the 
astrogeodetic  measurements  at  traverse  endpoints. 

The  Gauss-Markov  formula  relates  the  estimates  (de¬ 
flection  of  the  vertical  in  this  study)  to  the  separate  track 
data 


where 


lh\ 


-2 


\-N  / 


-  A  *-1 

=  0  d>  z 
TxzTzz  - 


(3.3-1) 
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$  =  cross-covariance  matrix  between  x  (truth 

quantities)  and  z 

d>  =  covariance  of  z 
zz  — 

x^  =  estimates  for  i*"*1  track 
z^  =  data  for  the  i1"*1  track 
N  =  total  number  of  tracks  involved. 

In  this  application  the  data  are  the  smoothed  estimates  of  the 
deflections  of  the  vertical  for  the  separate  tracks.  These 
quantities  are  "re-estimated"  based  on  data  from  nearby  tracks. 
The  covariance  of  the  estimation  errors  of  x  (x  =  x-x)  are 
determined  from 


A~~  —  <h 
YXX  XX 


)  d>  ^  0 
xz  Tzz  Tzx 


(3.3-2) 


where  $  is  the  a  priori  (reference)  covariance, 


To  describe  the  mechanics  of  the  estimation  procedure 
several  matrices  are  required.  The  matrix  6  relates  the 

XZ 

reference  field  deflections  to  the  single-track  deflection 
estimates 


I 

i 


♦xz  =  E<?5T) 


E<£lZ^)  £(£^2) .  •  -E^z^) 

£(*2-1^ 

...  .  E(x^Zjj) 


(3.3-3) 
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where  E  is  the  ensemble  expectation  operator.  The  matrix  4>zz 
is  the  covariance  matrix  of  the  single-track  deflection  estimates 


*zz  =  E<zzT)  = 


To  define  the  submatrices  needed  for  $  and  <t>zz »  expressions 
relating  the  reference  quantities  with  the  single-track  esti¬ 
mates  are  needed.  Two  forms  prove  to  be  useful  in  developing 
a  mathematically  and  computationally  tractable  formulation  of 
the  multiple  track  smoothing  problem.  The  first  is 

=  xi  +  g^  (Estimate  =  Truth  +  Error)  (3.3-5) 

where  g^  is  the  error  in  the  i1"*1  single-track  estimates. 
Through  optimal  smoothing  of  the  single-track  data  (including 
ZUPTS,  gradiometer  measurements,  or  star  sightings)  the  covar 
ia nee  between  the  estimate,  z^ ,  and  its  error,  g^ ,  is  zero. 
Thus,  the  covariance  equation  associated  with  Eq.  3.3-5  is 

Ri  =  I±  -  Pi  (3.3-6) 

where 

T  th 

R.  =  E(z.z.)  =  covariance  of  the  i  single-track 

1  11  estimates 

T 

I.  =  E(x.xT)  =  covariance  of  the  raw  gravity  field 
1  11  «. 

deflections  on  the  i  track 

rp  a.  L 

P.  =  E(g^.gT)  =  covariance  of  the  l  track  estimation 
1  11  errors 


E(£]£l)  E^z^) .  .  .E^zjlj) 


E<Z2-1 ^ 


E ) 


E(znzJ) 


(3.3-4) 
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A  state-space  covariance  simulation  of  the  errors  of  a  moving 
base  inertial  system  (a  fixed-point  smoothing  algorithm  is 
included)  is  used  to  generate  P^ .  The  reference  field  statis¬ 
tics,  1^,  are  known  a  priori  from  the  Markov  gravity  model. 

The  second  expression  (referred  to  as  the  "data  equa¬ 
tion")  relating  the  single-track  estimates  to  the  truth  values 
is 


z .  =  D.x .  +  e . 
-l  1-1  -l 


(3.3-7) 


where  e  is  an  error  that  is  orthogonal  to  x^  (the  covariance 
between  x^  and  e^  is  zero).  The  matrix  is  computed  from 


D.  =  I  -  P.lT1 
l  ii 


(3.3-8) 


where  I  is  the  identity  matrix. 

The  covariance  of  the  a  priori  deflections  may  be 
represented  by 


rxx 


1 

C12  • 

•  ‘  C1N 

21 

22  . 

r 

•  ’  zN 

N1 

.  .  . 

ZN 

(3.3-9) 


where  is  the  cross-covariance  of  the  deflections  between 
tracks  i  and  j .  The  cross  terms  approach  zero  as  the  tracks 
become  widely  separated. 


Using  the  expressions  listed  in  Eqs.  3.3-5  to  3.3-9, 
the  expectation  operations  in  Eqs.  3.3-3  and  3.3-4  are  carried 
out  to  yield 
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xz 


C 

C21U1 


cnidI 


C12  D2 


R, 


C1N°N 


VN 


(3.3-10) 


and 


zz 


R, 


D2C21D1 

DNCN1D1 


D1C12D2 


DC  D 
U1U1NUN 


R 


N 


(3.3-11) 


Note  that  since  the  between- track  correlations  of  the  a  priori 
deflections  approach  zero  as  the  track  separation  increases, 
the  dimensionality  of  the  problem  can  be  maintained  at  a  tract 
able  level  by  incorporating  data  from  nearby  tracks  only. 
Processing  data  from  remote  tracks  offers  little  accuracy 
improvement.  Note  that  appropriate  quantification  of  "remote" 
can  be  established  by  simulation.  Examples  are  provided  in 
the  next  section. 


3. A  NUMERICAL  RESULTS 

The  improvement  in  deflection  estimation  accuracy 
achievable  by  combining  single-track  deflection  of  the  vertical 
estimates  from  moving  base  inertial  systems,  along  with  the 
astro-measurements,  is  reported  here.  Such  results  are  indi¬ 
cative  of  the  benefits  derived  from  a  post  survey  regional 
adjustment  (data  reduction  process)  based  on  the  formulation 
described  in  the  previous  section.  The  example  case  used  in¬ 
volves  two  parallel  tracks  with  the  separation  distance  varied 
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Figure  3.4-1  Accuracy  Potential  of  Conventional  Inertial 

System  Two-Track  Survey 

as  a  parameter.  Accuracy  results  for  the  three  systems  simu¬ 
lated,  a  conventional  inertial  system,  a  gradiometer-aided 
system,  and  the  NAS-26  astroinertial  system,  are  presented  and 
discussed.  Deflection  of  the  vertical  estimation  accuracy  is 
presented  in  graphs  which  also  indicate  the  corresponding  single¬ 
traverse  accuracy  found  in  Ref.  1.  North  deflection  of  the 
vertical  post-survey  residual  rms  error  performance  (the  direc¬ 
tion  transverse  to  the  track  direction)  is  presented  in  the 
following  paragraph.  Improvement  in  the  along- track  components 
of  the  deflection  (in  this  case,  East)  with  multiple  survey 
tracks,  is  somewhat  smaller. 

Note,  in  Fig.  3.4-1,  that  no  improvement  occurs  at 
the  traverse  endpoints.  This  is  to  be  expected  for  the  single- 
track  traverse  since  the  inertial  survey  system  provides  relative 
deflection  information,  i.e.,  the  deflection  change  from  an 
endpoint  value.  Although  a  completely  error  free  survey  would 
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link  endpoints  as  two  effectively  collocated  measurements,  the 
inertial  system  accuracies  corresponding  to  the  results  in 
Fig.  3.4-1  are  insufficient  to  provide  a  significant  amount  of 
endpoint  coupling.  A  counter  example,  illustrating  a  case 
wherein  the  effective  endpoint  coupling  is  noticeable,  is  pre¬ 
sented  in  Fig.  3.4-2. 

For  the  multiple  track  cases  in  Fig.  3.4-1,  the  0.2  sec 
endpoint  deflection  accuracies  also  do  not  improve.  This  is 
due  to  the  recognition  that  a  vertical  deflection  station  of 
this  accuracy  would  usually  be  the  result  of  a  previous  regional 
adjustment  of  deflection  data.  This  was  taken  into  account  in 
the  simulation  study  by  treating  the  endpoint  deflection  errors 
for  adjacent  tracks  as  fully  correlated. 

For  the  estimation  accuracy  potential  of  a  conventional 
moving  base  inertial  system  presented  in  Fig.  3.4-1,  the  base¬ 
line  gyro  and  accelerometer  error  models,  summarized  in  Table 
3.2-3,  are  used.  ZUPTs  are  performed  at  4  minute  intervals 
until  the  end  of  the  survey.  Figure  3.4-1  indicates  the  de¬ 
flection  of  the  vertical  estimation  uncertainty  as  a  result  of 
optimally  combining  data  from  two  parallel  tracks.  This  example 
is  chosen  for  its  simplicity  and  illustrates  the  variation  of 
accuracy  with  track  separation.  The  processing  method,  however, 
is  applicable  for  any  reasonable  number  of  tracks  with  no  inher¬ 
ent  geometry  restraints.  The  tracks  run  from  West  to  East  so 
that  the  North  deflection  is  the  cross-track  component.  Note 
that  the  smoothing  operation,  which  utilizes  all  survey  measure¬ 
ments  in  forming  an  estimate  at  any  point  in  the  survey,  ensures 
a  near  symmetric  error  profile.  The  highest  uncertainty  occurs 
at  the  middle  of  the  traverse,  farthest  from  the  calibrated 
endpoints.  The  results  show  a  significant  improvement  with 
respect  to  the  single-traverse  accuracy  depending  on  the  track 
separation  distance.  Closer  track  spacing  yields  better  accu¬ 
racy.  Also,  the  degradation  in  accuracy  as  the  tracks  are 
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separated  is  clearly  seen.  This  is  due  to  the  decreasing  spatial 
correlation  of  the  gravity  field. 

The  Bell  gradiometer  error  model  is  integrated  with 
the  baseline  navigation  system  and  sensor  error  models  to  simu¬ 
late  an  astrogeodetic/gradiometer-aided  inertial  system  survey 
scenario.  Endpoint  calibration  data,  survey  length,  cruise 
velocity,  and  the  gravity  disturbance  model  remain  unchanged. 

No  ZUPTs  are  performed.  Note  from  the  performance  indicated 
in  Fig.  3.4-2  that  ZUPTs  are  not  needed  for  deflection  survey- 
related  gradiometer  operations. 
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Figure  3.4-2  Accuracy  Potential  of  Gradiometer-Aided 

Inertial  System  Two-Track  Survey 
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In  computing  the  post  survey  deflection  of  the  ver¬ 
tical  error  residuals  for  the  parallel  track  geometry,  it  was 
found  that  the  single- traverse  accuracy  derived  from  the  gra- 
diometer-aided  system  could  not  be  significantly  improved  by 
mixing  data  from  tracks  as  close  as  8.95  km,  one-half  the  grav¬ 
ity  disturbance  characteristic  distance.  This  is  shown  in 
Fig.  3.4-2. 

The  gradiometer,  being  a  highly  accurate  and  precise 
interpolator,  in  effect  carries  the  astrogeodetic  endpoint 
measurements  along  the  entire  traverse  with  a  minimum  of  error. 
Little  is  gained  by  performing  post  survey  regional  adjustments 
on  data  collected  by  gradiometer-aided  inertial  systems.  Of 
course,  adjustment  of  the  endpoint  deflection  stations  according 
to  the  procedures  outlined  in  Ref.  2  could  offer  the  additional 
improvements  indicated  by  a  separate  regional  deflection  station 
adjustment . 

The  survey  scenario  for  the  NAS-26  simulation  mimicked 
that  for  the  conventional  inertial  system  case.  ZUPTs  were  per¬ 
formed  at  4  minute  intervals.  However,  star  sightings  were 
also  modeled  and  took  place  every  20  seconds.  The  accuracy 
improvement  obtained  by  combining  deflection  estimates  is  pre¬ 
sented  in  Fig.  3.4-3.  In  comparison  with  the  conventional 
system  results  of  Fig.  3.4-1  the  NAS-26  single-track  estimates 
are  better.  This  is  a  result  of  the  ability  of  the  NAS-26  to 
finely  calibrate  the  gyros  in  real  time  using  the  star  sighting 
information.  Since  gyro  errors  are  a  major  contributor  to  the 
bowed  error  profile  of  Fig.  3.4-1,  as  reported  in  Ref.  1,  the 
elimination  of  such  errors  results  in  the  flattened  curve  in 
Fig.  3.4-3.  Accordingly,  the  accuracy  improvement  realized  by 
combining  NAS-26  data  from  two  parallel  tracks  is  not  as  marked 
as  that  observed  from  the  conventional  inertial  system.  This 
is  not  surprising  since  there  is  less  error  to  be  improved 
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upon.  Using  the  error  models  and  survey  scenario  described  in 
this  effort,  the  NAS-26  astroinertial  navigator  is  found  to 
yield  deflection  estimates  of  intermediate  accuracy  in  compari¬ 
son  to  that  obtained  by  conventional  and  gradiometer-aided 
inertial  systems. 
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4. 


COSTS  AND  LOGISTICS 


In  this  chapter  assessments  of  relative  costs  and 
operational  considerations,  such  as  survey  time,  are  presented 
for  the  system  configurations  analyzed.  Accuracy  comparisons 
are  made  using  the  results  presented  in  previous  chapters. 

Astrogeodetic  measurements  of  the  deflection  of  the 
vertical  must  be  made  to  achieve  the  high  accuracy  (0.2  sec) 
desired  from  modern  surveys.  The  measurements  are  time- 
consuming,  however,  requiring  a  skilled  crew  to  take  one  night 
or  more  per  measurement.  Therefore,  such  data  is  intended  to 
supplement  other  measurement  procedures  and  aid  in  lowering 
deflection  of  the  vertical  residual  errors.  In  Chapter  2  it 
was  shown  that  gravimeter  measurements,  which  are  simpler  and 
less  time-consuming  to  make,  combined  with  astrogeodetic  meas¬ 
urements,  yield  estimation  uncertainties  dependent  on  the 
spacing  of  the  astrogeodetic  measurements  and  the  number  of 
gravimeter  measurements,  with  the  potential  of  achieving  high 
accuracies.  Optimal  processing  methodology  is  used  in  assess¬ 
ing  the  tradeoffs.  Since  astrogeodetic  and  gravimetry  tech¬ 
nology  is  mature,  little  technical  risk  is  incurred  by  taking 
this  approach  to  survey  a  given  area.  Incremental  productivity 
gains  may  be  realized. 

Inertial  and  gradiometric  surveys  performed  on  the 
ground  or  in  the  air  offer  higher  survey  productivity  and  can 
utilize  infrequent  astrogeodetic  "benchmark”  data  in  the  process 
ing  algorithms.  Despite  optimal  processing  and  post  survey 
regional  adjustments  involving  two  tracks,  conventional  inertial 
surveys  have  difficulty  reaching  deflection  of  the  vertical 
estimation  accuracies  of  0.2  sec.  Inertial  systems,  however, 
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have  reached  a  high  level  of  positioning  precision  and  are 
readily  available.  As  mentioned  in  Ref.  1  such  systems  could 
be  used  to  provide  accurate  positional  data  for  surveys  involv¬ 
ing  a  gradiometer,  which  can  provide  very  accurate  deflection 
of  the  vertical  estimates  given  astrogeodetic  endpoint  data. 

High  development  and  procurement  costs  are  associated  with 
gradiometers ,  and  they  have  yet  to  be  fully  demonstrated  in 
prototype  testing.  Therefore,  greater  technical  risk  follows 
the  use  of  a  gradiometer. 

The  Northrop  NAS-26  astroinertial  navigation  system, 
presently  not  used  for  other  than  experimental  gravity  data 
surveying,  appears  to  estimate  deflections  of  the  vertical 
well.  Marginal  improvements  in  estimation  accuracy  are  found 
when  analyzing  a  two- track  survey.  The  successful  operation 
of  the  NAS-26  is  dependent  on  weather  conditions  (clear  skies), 
and  added  software  would  be  necessary  to  handle  ZUPTs  and  gravity 
data  processing. 

Figure  4-1  sketches  the  relative  costs  vs  accuracy  of 
various  multisensor  combinations.  Multiple  tracks  are  illus¬ 
trated  to  increase  survey  costs  for  a  given  system  configura¬ 
tion  (through  increased  survey  effect).  The  gradiometer  con¬ 
figuration  exhibits  the  highest  accuracy,  with  correspondingly 
high  costs.  The  astrogravimetric  approach  could  be  especially 
attractive  in  areas  where  considerable  anomaly  data  is  available. 
Further  study  identifying  specific  candidate  locations  for 
test  surveys  or  astrogravimetric  processing  of  existing  data 
is  needed. 

The  trends  illustrated  in  Fig.  4-1  provide  indications 
of  the  relative  costs  of  different  options  for  precise  deflec¬ 
tion  of  the  vertical  surveying.  Detailed  consideration  of 
cost  and  logistical  factors,  of  course,  is  a  multidimensional 
problem  in  its  own  right,  with  few  overall  value  measures  which 
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are  universally  accepted.  Although  a  complete  and  detailed 
study  of  all  such  factors  is  beyond  the  resources  of  this  study 
Tables  4-1  and  4-2  present  several  significant  considerations 
for  each  class  of  surveying  approach.  Of  course  each  option 
listed  in  the  table  could  be  pursued  in  several  variant  forms. 
For  example,  the  different  remedial  procedures  and  processes 
discussed  in  Ref.  2  can  be  applied  to  T-4  Theodolite-based 
astrogeodetic  surveys.  Similarly,  inertial  surveying  system 
gyros  and  accelerometers  can  be  replaced  with  more  sensitive 
components.  While  such  refinements  of  the  main  technological 
thrusts  in  deflection  surveying  are  important,  they  typically 
provide  single  benefits  such  as  more  accuracy,  lower  cost,  more 
efficiency,  etc.  Since  the  incremental  improvement  offered  by 
such  refinement  technology  is  usually  straightforward  to  assess 
Tables  4-1  and  4-2  instead  address  alternative  approaches  which 
span  different  technologies. 
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DATA 

TYPE  (.AFTER 


APPROXIMATE 
SYSTEM 


SURVEY  APPROACHES  REDUCTION)  HARDWARE  INVOLVED  (sec) 
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TABLE  4-1 

COST  AND  LOGISTICS  CONSIDERATIONS  FOR 
ALTERNATIVE  FIXED  STATION  TECHNIQUES 
TO  DETERMINE  DEFLECTIONS 
OF  THE  VERTICAL 
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TABLE  4-2 

COST  AND  LOGISTICS  CONSIDERATIONS  FOR 
ALTERNATIVE  MOVING  BASE  TECHNIQUES 
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4.1  CONVENTIONAL  ASTROGEODETIC  SURVEY 

Conventional  theodolite  surveys,  discussed  in  detail 
in  Ref.  6,  are  the  current  mainstay  of  deflection  of  the  ver¬ 
tical  determination  in  the  United  States.  The  technology  is 
mature  and  the  tradeoffs  among  alternative  data  reduction 
methodologies  are  well  understood.  The  key  disadvantage  with 
the  astrogeodetic  method  is  the  inherent  difficulty  of  using  a 
theodolite  to  provide  astronometric  data  consistent  with  0.2 
arcsecond  accuracy.  This  element  is  indicated  in  Table  4-1 
both  in  the  time  required  to  train  surveying  personnel  (note 
that  many  in  the  community  view  several  months  as  an  optimistic 
estimate  of  the  apprenticeship  period  for  a  first-order  astro 
surveyor)  and  the  time  and  cost  needed  to  obtain  a  single  astro¬ 
geodetic  station.  Cost  estimates  given  in  the  table  are  based 
on  discussions  with  members  of  the  surveying  community. 


4.2  ASTROLABE -BAS ED  DEFLECTION  SURVEYS 

Recent  interest  in  using  astrolabes  for  deflection 
determination  has  arisen  primarily  as  a  response  to  the  train¬ 
ing  and  observation  time  requirements  for  conventional  theo¬ 
dolite  surveys.  The  key  justification  for  the  consideration 
of  astrolabes  is  the  more  relaxed  training  and  observing  regime. 
In  Table  4-1,  a  Danjon  astrolabe  is  considered  as  typifying 
field  portable  instrumentation.  "Field  portable"  is  used  in 
the  same  sense  as  often  applied  to  the  Wild  T-4  theodolite. 
Although  not  easily  handled,  a  party  of  two  can  manage  the 
equipment  on  all  except  the  most  inhospitable  terrain.  Pricing 
information  is  estimated  with  the  note  that  the  Danjon  astrolobe 
is  no  longer  made,  and  no  predictable  second-hand  market  exists. 
However  discussions  with  various  user  community  members  indicate 
the  amount  given  in  Table  4-1  is  a  reasonable  guess  at  costs. 
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Since  the  number  of  astrolabes  available  at  any  given 
time  is  likely  to  be  low,  any  large  scale  or  quick-response 
survey  requirement  which  would  be  addressed  using  astrolabes 
would  likely  be  subject  to  manufacturing  or  procurement  delays 
of  several  years.  Alternatively  if  custom  instruments  were 
built,  or  a  new  production  facility  opened,  the  cost  of  a  dozen 
or  so  astrolabes  produced  on  an  accelerated  basis  could  easily 
reach  a  unit  cost  more  than  an  order  of  magnitude  higher  than 
given  in  Table  4-1. 

4.3  T-4  WITH  CHARGE  COUPLED  DEVICE  (CCD)  EYEPIECE 

The  CCD  eyepice  is  a  technology  which  promises  modest 
accuracy  improvements  (reduces  human  observer-induced  measure¬ 
ment  uncertainty)  and  significant  productivity  increases. 

Costs  indicated  in  Table  4-1  would  add  to  amounts  needed  for 
the  theodolite  and  yet-to-be  developed  production  survey  computer 
software.  The  decreased  cost  per  data  point  (versus  conventional 
astrosurvey)  indicated  in  Table  4-1,  is  estimated  to  be  gained 
largely  through  data  processing  efficiencies  and  reduction  in 
the  amount  of  rejected  data. 


4.4  ASTROLABE  WITH  CCD  EYEPIECE  AND  TWO  COLOR 
REFRACTOMETER  (TCR) 

This  developmental  technology  promises  significant 
accuracy  and  productivity  improvements  beyond  current  T-4  per¬ 
formance  through  the  elimination  of  observer  uncertainty,  con¬ 
trol  of  refraction  effects  and  utilization  of  the  astrolabe's 
efficiency.  The  system  cost  estimates  of  Table  4-1  are  based 
on  the  procurement  of  hardware  for  a  prototype;  they  exclude 
costs  for  development  of  supporting  computer  software  and  test¬ 
ing  costs. 
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4.5  GRAVIMETRIC  SURVEY 

Gravimetric  surveys,  while  not  generally  as  accurate 
for  deflection  of  the  vertical  generation  as  astrogeodetic 
techniques,  are  attractive  because  of  the  relatively  low  equip¬ 
ment  costs,  ease  of  use  and  survey  productivity  limited  essen¬ 
tially  by  how  fast  the  gravimeter  can  be  transported  from  site 
to  site.  Of  course,  although  the  cost  per  gravity  anomaly 
point  is  low,  it  must  be  kept  in  mind  that  the  processing  effort 
required  to  transform  gravity  data  into  deflections  is  consider¬ 
able.  For  subarcsecond  accuracy,  not  only  are  many  gravimetric 
measurements  required  near  each  deflection  point  (inner  zone 
densification)  but  the  processing  must  be  coordinated  with 
regional  and  worldwide  data  base  holdings  to  properly  account 
for  distant  zone  weighting.  Although  exceptions  have  been 
noted  (see  Ref.  2)  ,  gravimetry  alone  is  generally  not  con¬ 
sidered  viable  for  high  accuracy  deflection  determination  in 
the  CONUS. 

4.6  ASTROGRAVIMETRIC  SURVEY 

The  key  promise  of  the  astrogravimetric  method  is  the 
relaxation  of  astro  densification  requirements  through  inter¬ 
polation  using  less  costly  gravimetric  measurements.  As  indi¬ 
cated  in  Table  4-1,  the  results  of  the  present  study  indicate 
that  desired  area-wide  survey  accuracies  can  be  achieved  with 
an  astrogravimetric  approach.  However,  less  clear  is  whether, 
in  actual  survey  practice,  astrogravimetric  surveying  would 
offer  a  significant  cost  or  time  savings  over  an  astro  survey 
of  equivalent  accuracy.  Definitive  answers  could  be  provided 


Also,  currently-used  launch  region  gravity  model  techniques 
for  certain  land-based  ICBMs  utilize  a  variant  of  the  gravi 
metric  method. 
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either  by  an  intensive  directed  study  or  implementation  of  a 
proto  survey.  The  latter  would  require  a  significant  commit¬ 
ment  to  computer  software  development  for  data  integration. 
The  indicated  acquisition  costs  assume  the  use  of  one  gravi¬ 
meter  and  one  T-4  or  equivalent  theodolite. 


4.7  SURFACE  OR  HELICOPTER  INERTIAL  SURVEY 

This  category  in  Table  4-2  describes  current  inertial 
survey  technology.  Costs  per  track  of  data  are  driven  prima¬ 
rily  by  field  crew  and  equipment  maintenance  costs,  (assuming 
high  annual  utilization  rates).  The  key  advantage  of  inertial 
survey  approaches,  or  any  other  moving-base  technique,  is  rapid 
data  gathering,  amenability  to  automation  and  a  minimal  need 
for  personnel  with  highly  specialized  skills.  Note  that,  while 
the  accuracy  indicated  in  Table  4-2  does  not  meet  the  0.2  sec 
goal ,  additional  improvements  to  inertial  survey  systems  are 
possible.  These  are  detailed  and  evaluated  in  Ref.  1.  One 
such  improvement  is  discussed  in  the  next  paragraph. 


4.8  INERTIAL  SURVEY  USING  A  STAR  TRACKER 


The  only  currently-available  system  utilizing  this 
aided-inertial  configuration  is  the  NAS-26  prototype.  Although 
some  development  would  be  required  to  tranform  this  system 
into  a  production  survey  configuration,  certain  test  results 
as  well  as  the  simulation  studies  described  earlier  in  this 
report  are  encouraging.  As  indicated  in  Table  4-2,  current 


For  example,  taking  a  five  year  amortization  of  a  one  million 
dollar  acquisition  and  spares  cost,  adding  fifty  thousand 
dollars  of  average  annual  finance  charges  and  assuming  100  traverses 
per  year  of  utilization,  gives  a  capital  cost  alone  of  $250 
per  traverse. 
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hardware  costs  would  decrease  once  the  system,  now  baselined 
for  B-l  test  instrumentation,  is  in  full  production.  Although 
the  survey  accuracy  indicated  in  Table  4-2  is  not  at  goal  levels, 
improvements  such  as  those  discussed  in  Ref.  1  are  possible. 

One  particularly  effective  expedient,  of  course,  is  to  increase, 
to  the  extent  possible,  deflection  calibration  accuracy  of 
traverse  terminal  points. 


4.9  INERTIAL/GRAD IOMETER  SURVEY 

The  gravity  gradiometer,  now  under  development  with 
DMA  sponsorship,  has  already  demonstrated  laboratory  performance 
sufficient  to  reach  the  accuracy  levels  indicated  in  Table  4-1. 
At-sea  tests  have  verified  ruggedness.  However,  several  years 
of  applied  development  await  production  survey  utilization  of 
gradiometer  systems.  Nonetheless,  when  operational,  the  econom¬ 
ics  and  logistics  of  their  utilization  should  be  nearly  equiva¬ 
lent  to  those  of  conventional  inertial  survey  equipment. 


4.10  FIXED  WING  AIRCRAFT  GRADIOMETER  SURVEY 

The  long  term  potential  of  the  gradiometer  as  a  rapid 
survey  tool  is  efficiently  exploited  by  using  a  fixed  wing 
aircraft  as  the  host  vehicle.  Although  well  into  the  future, 
current  data  suggest  that,  for  surveying  very  large  regions 
quickly,  the  very  high  development  and  acquisition  costs  for 
an  airborne  gradiometric  survey  system  could  be  worthwhile. 

Key  technical  risks  lie  in  the  areas  of  successful  system  in¬ 
tegration  of  existing  technology  and  aircraft,  and  the  imple¬ 
mentation  of  techniques  to  process  the  very  large  volumes  of 
data  returned  by  an  airborne  gradiometric  sensor.  An  addi¬ 
tional  potential  benefit  of  an  airborne  gradiometer  system  is 
its  applicability  to  broad  ocean  areas  as  well  as  over  land. 
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4.11  OVERALL  PERSPECTIVE 


Tables  4-1  and  4-2  supports  several  observations 


•  Technologically  advanced  alternatives  to 
conventional  astrogeodetic  surveys  can 
provide  considerably  increased  deflection 
data  acquisition  capaoility  at  current 
first-order  accuracies  or  better.  How¬ 
ever,  these  alternatives  have  associated 
research  and  acquisition  costs  which 
tend  to  increase  exponentially  with  capa¬ 
bility.  These  costs  may  be  justifiable 
in  terms  of  per  unit  deflection  data 
cost  if  alternative  new  systems  are 
heavily  utilized. 

•  Astrogeodetic  survey  techniques,  especially 
those  involving  astrolabe  use,  updated 

to  include  automated  observing  and  more 
computerized  data  reduction  offer  marked 
relief  from  traditional  difficulties  in 
training  and  data  quality  assurance. 

Order  of  magnitude  improvements  in  data 
acquisition  rates,  however,  are  unlikely 
due  to  star  viewing  time  constraints. 

•  The  costs  of  developing  and  maintaining 
an  automated  or  semiautomated  astrogeo¬ 
detic  technology  base  are  modest  compared 
with  those  associated  with  advancing 
moving-base  survey  techniques. 

•  Economic  benefits  associated  with  util¬ 
izing  astrogravimetric  methods  may  be 
realizable  but  modest. 

•  Moving-base  methods,  the  only  feasible 
approaches  to  densifying  deflection  data 
over  large  areas  (e.g.,  all  CONUS  coast¬ 
lines,  new  multi-state  areas,  etc.)  within 
a  reasonable  timespan,  are  within  reach 
using  presently-available  technology. 

The  cost  and  risk  of  this  technique  vary 
in  approximate  proportion  to  the  ultimate 
survey  accuracy  and  speed  offered. 
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•  The  recurrent  technical  risk  theme  in 
most  alternative  deflection  measurement 
schemes  is  the  development  of  computer 
software.  Although  the  software  needs  of 
most  candiate  sensors  are  easily  addressed, 
little  has  been  done  toward  integrating 

the  widespread  survey  data  (or  even  organiz¬ 
ing  its  acquisition)  which  will  be  provided 
by  new  alternative  techniques.  These 
new  measurements,  which  will  be  obtained 
on  a  large  scale  at  sub  arcsecond  accuracies, 
must  be  integrated  into  existing  data 
bases  without  loss  of  information.  In 
addition ,  this  processing  must  be  accom¬ 
plished  without  overloading  the  mapping 
community's  ability  to  perform  the  neces¬ 
sary  data  reductions.  Software  and  data 
base  techniques,  similar  to  those  embodied 
in  the  present  DoD  gravity  anomaly  library 
and  DMA  weapon  support  systems,  are  needed 
at  the  instrument/field  level  to  prepare 
for  the  future  availability  of  very  accu¬ 
rate  and  highly  densified  deflection  of 
the  vertical  measurements. 

•  A  relatively  low  cost  (when  compared  to 
the  cost  of  moving  base  survey  hardware 
development)  interim  approach  to  main¬ 
taining  a  deflection  survey  "breakout" 
capability  for  the  next  five-to-ten  years 
would  be  to  simply  proliferate  astrolabe- 
trained  observers,  acquire  instruments 
and  initiate  a  data  acquisition  program. 

Such  an  approach  could  adequately  address 
deflection  of  the  vertical  requirements 
for  strategic  weapons  systems  where  the 
data  needs  are  restricted  to  limited 
areas  and  also  for  many  tactical  needs. 
Examples  include  deflection  initializa¬ 
tion  stations,  aircraft  inertial  systems, 
vertical  deflection  references  at  govern¬ 
ment  test  sites  and  independent  deflec¬ 
tion  sites  used  as  survey  tie  points. 
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5. 


SUMMARY  AND  CONCLUSIONS 


Alternative  multisensor  deflection  of  the  vertical 
survey  approaches  were  examined  which  offer  promise  of  ap¬ 
proaching  or  bettering  an  accuracy  goal  of  0.2  sec.  Specific 
techniques  studied  and  their  projected  performance  include: 


•  Gravimetric  Enhancement  of  Astrogeodetic 
Data  Fields  -  The  astrogravimetric  approach 
offers  a  survey  design  "tradeoff"  in 
which  many  gravimetric  measurements  can 

be  substituted  for  a  few  astrogeodetic 
stations.  However,  a  minimum  number  of 
astro  measurements  is  always  required  to 
reach  deflection  accuracy  goal  levels. 

This  number  varies  with  proximity  of  the 
deflection  estimate  to  a  location  where 
a  deflection  observation  is  available. 

The  astro  site  spacing  also  varies  with 
the  accuracy  of  both  the  astrogeodetic 
and  gravimetric  data  and  the  correlation 
distance  of  the  gravity  field.  Typical 
deflection  and  anomaly  measurement  densi¬ 
ties,  which  offer  goal  level  deflection 
accuracy,  range  from  no  gravity  anomaly 
measurements  and  an  astro  site  spacing 
of  5  km  to  an  anomaly  spacing  of  1.7  km 
with  an  astro  spacing  of  15  km. 

•  Regional  Adjustment  of  Conventional  Iner- 
Fial  Survey  Traverses  Using  Astrogeodetic 
Endpoint  Data  -  Post-survey  adjustment 

of  inertial  survey  measurements,  using 
optimum  smoothing  techniques,  offers 
deflection  accuracy  improvements  over 
single  track  traverses  up  to  20%  per 
additional  data  track.  However  the  ac¬ 
complishment  of  goal  level  accuracies  is 
still  difficult.  Using  conventional 
inertial  technology  (e.g.,  RGSS)  and  two 
35  km  parallel  tracks,  mid-traverse  de¬ 
flection  recovery  is  0.6  sec  rms. 
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•  Regional  Adjustment  of  NAS-26  Quality, 
Tnertial  Survey  Traverses  -  Post-survey 
adjustment  of  NAS-26  survey  data  offers 
modest  improvement  beyond  single  track 
accuracies.  For  traverses  between  15 
and  35  km  long,  typical  midpoint  deflec¬ 
tion  accuracy  is  near  0.4  sec  rms, 

•  Regional  adjustment  of  Gravity  Gradio¬ 
meter  Survey  Traverses  -  The  gradiometer 
effectively  replicates  the  deflection 
endpoint  accuracy  at  all  traverse  points. 
Regional  adjustment  of  gradiometer  trav¬ 
erse  data  provides  negligible  deflection 
recovery  improvement  behond  that  obtain¬ 
able  by  optimally  smoothing  each  traverse 


The  formulation  for  combining  multisensor/multitrack 
data  from  moving  base  inertial  system  surveys  for  gravity  data 
estimation  is  a  major  product  of  this  effort.  The  methodology 
is  general  in  terms  of  allowable  measurement  types  and  track 
arrangements.  Future  work  should  focus  on  such  cases  as  cross¬ 
ing  track  and  repeat  traverse  scenarios.  Most  importantly, 
the  methodology  should  be  used  to  process  actual  RGSS  data  ro 
realize,  at  the  earliest  possible  date,  the  benefits  of  the 
approach . 


The  information  presented  in  Tables  4-1  and  4-2  pro¬ 
vides  a  capsule  summary  of  major  deflection  survey  alternatives. 
These  alternatives  include  both  single  optical  sensor  approaches 
as  well  as  integrated  multisensor  methods. 

This  study,  as  well  as  two  earlier  related  studies 
(Refs.  1  and  2),  supports  the  following  conclusions: 

•  Technologically  advanced  deflection  of 
the  vertical  survey  systems  now  under 
development  offer  promise  of  0.2  sec 
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accuracy  or  better,  as  well  as  relief 
from  the  time  constraints  now  associated 
with  acquiring  deflection  data. 

•  Only  modest  relief  can  be  expected  from 
the  current  cost  of  high  accuracy  deflec¬ 
tion  of  the  vertical  stations. 

•  Efficient  utilization  of  new  deflection 
of  the  vertical  survey  approaches  accents 
the  need  to  direct  resources  toward  auto¬ 
mation.  Particular  attention  to  the  devel¬ 
opment  and  operation  of  software  and 
computer  support  systems  will  be  especially 
needed . 


Further  study  is  recommended  toward  the  effective 
organization  and  utilization  of  deflection  of  the  vertical 
data.  Considerable  new  amounts  of  data  of  sub-arcsecond 
accuracy  will  become  available  as  the  techniques  described 
in  this  and  earlier  reports  (e.g.,  Refs.  1  and  2)  are  im¬ 
plemented  . 
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APPENDIX  A 

SEQUENTIAL  COLLOCATION  FORMULAS 


Sequential  collocation  offers  a  convenient  method  of 
processing  measurement  data  from  one  sensor  at  a  time  in  an 
estimation  problem  in  which  several  sensors  may  be  involved. 
Processing  multisensor  data  sequentially  reveals  the  impact 
each  sensor  has  on  estimation  accuracy.  Complex  iteration 
algorithms  are  available  in  the  literature  (Ref.  5)  for  imple¬ 
menting  the  sequential  collocation  technique.  For  convenience 
the  sequential  collocation  estimation  and  estimation  error 
covariance  formulas  are  listed  in  this  Appendix  for  three  meas¬ 
urement  types. 


A. 1  PROBLEM  FORMULATION 

It  is  desired  to  estimate  a  vector  quantity  X  given 
measurements  of  u,  v,  and  w.  A  measurement  of  X  itself  may  be 
made,  and  correlations  may,  of  course,  exist  among  any  of  the 
quantities.  The  dimensions  of  the  vectors  need  not  be  equal. 

As  an  example,  suppose  it  is  required  to  estimate  deflections 
of  the  vertical  in  a  given  survey  area  wherein  gravity  anomaly 
and  undulation  measurements  are  available.  These  three  quanti¬ 
ties  are  correlated  to  some  degree,  and  an  appropriate  gravity 
disturbance  model  is  used  to  describe  the  covariance  structure 
of  the  variables.  For  example,  collocation  formulas  could  be 
used  to  estimate  the  deflections  and  their  associated  uncer¬ 
tainty  using  the  anomaly  measurements  alone  and  then  in  con¬ 
junction  with  undulation  measurements. 
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Table  A-l  defines  the  notation  used  in  the  formula- 
T 

tion.  Note  that  C  ,  =  C,  . 

ab  ba 

TABLE  A-l 

NOTATION  AND  NOMENCLATURE  FOR 
COLLOCATION  FORMULAS 


-abc  =  est:*-mate  °f  X  given  measurements  of  a,  b,  and  c 

Za  =  noisy  measurement  of  a  (a  +  rj.  ) 

R  =  covariance  of  n 
a  -La 

Cab  C(je  =  covariance  of  a  and  b  given  measurements  of  c,  d, 
13  and  e . 


A. 2  SEQUENTIAL  COLLOCATION  FORMULAS 

In  Table  A-2  the  formulas  for  the  first  estimate  of  X 
and  its  associated  error  covariance  matrix,  given  the  measure¬ 
ments  Zu  of  the  first  sensor,  are  presented.  These  are  famil¬ 
iar  expressions  in  minimum  variance  estimation.  Given  the 
measurements  Zv  of  the  second  sensor,  the  estimate  and  its 
error  covariance  matrix  are  calculated  as  shown  in  Table  A-3. 
Note  that  in  a  state-space  format  new  measurement  (H)  and 
measurement  noise  (R)  matrices  would  be  brought  in  to  perform 
the  update.  However,  using  collocation,  the  computation  of 
conditional  covariance  matrices  is  necessary.  In  particular, 
the  results  of  the  first  update  are  needed  in  the  computation 
of  the  second.  Continuing  the  process,  the  results  of  util¬ 
izing  Zv  are  embodied  in  the  formulas  of  Table  A-4.  These 
formulas  process  the  final  measurements,  Z  . 
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TABLE  A- 2 

COLLOCATION  FORMULAS  FOR 
FIRST  MEASUREMENT  TYPE 


X 

=  C 

1C  +  R  ]" 

uu  u 

1 

Z 

— u 

XU 

— u 

C 

=  C 

-  C  [C 

+ 

R  l'1  C 
u  ux 

XX  .  u 

XX 

XU  uu 

TABLE  A- 3 

COLLOCATION  FORMULAS  FOR 
SECOND  MEASUREMENT  TYPE 


=  X  +  C  (C  +  R  )_1IZ  -C  [C  +  R  )~lZ  ] 

— u  xv . u  vv . u  v  -v  vu 1  uu  u  J  -u 1 


xx .  uv 


Cxx.u  ”  Cxv.u^Cw.u  +  V  Cvx.u 


where 


XV  .  u 

=  C 

-  C 

[c 

+  R  l"1  C 

U 1  UV 

XV 

XU 

UU 

'w.  u 

s  c  , 

-  c 

[C 

+  R  J’1  C 

U  UV 

vv 

VU 

1  UU 

TABLE  A- 4 

COLLOCATION  FORMULAS  FOR 
THIRD  MEASUREMENT  TYPE 


A  A 

X  =  X 

+  c 

(C  +  R 

)_1  Z  -C  (C  +  R  )~1[z  -C  (C  +  R  )~lz  ] 

— UVW  “UV  xv.uv  wv. uv  w 

—w  WV.U  W.U  V  V  VU  uu  u  — u 

c  =  c 

C  [C  + 

R  J_1  C 

XX  .  UVW  XX  .  uv 

XW .  uv  ww .  uv 

w  wx . uv 

where, 

C  = 

c 

-  C  (C  + 

R  )-1  C 

XW.  uv 

xw.u 

xv .  u  w .  u 

V  vw.  u 

C  = 

C 

-  C  (C  + 

R  )_1  C 

WW.  uv 

WW.U 

wv.u  w.u 

v  vw.u 

c  = 

c 

C  (C  +  R  )_1 

C 

xw.u 

XW 

xu  uu  u 

uw 

C  = 

C 

c  (c  +  r  r1 

C 

WV.  u 

WV 

wu  uu  u 

uv 

c  = 

c 

C  (C  +  R  )_1 

C 

ww.u 

WW 

wu  uu  u 

UW  j 
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The  above  sequential  collocation  procedure  is  compu¬ 
tationally  equivalent  to  the  straightforward  batch  formulation 
of  Table  A-5.  In  Table  A-5  all  sensor  measurements  are  proc¬ 
essed  at  once  to  yield  the  final  estimate  vector  and  its  error 
covariance  matrix.  The  effects  of  the  individual  sensors  can¬ 
not  be  seen.  Although  a  large  matrix  must  be  inverted  in  the 
batch  scheme  (and  precludes  using  this  approach  for  a  very 
large  number  of  measurements) '  no  real  computational  savings 
are  realized  in  the  sequential  algorithm  due  to  the  added  com¬ 
plexity  of  computing  the  conditional  covariance  matrices.  Its 
main  advantage  is  the  conceptual,  sensor-specific  clarity 
it  offers. 


TABLE  A-5 

COLLOCATION  BATCH  SCHEME  FOR 
THREE  MEASUREMENT  TYPES 


X  = 

(c 

C  C  ] 

C 

+R  C  C  _1  Z 

-uvw 

1  XU 

XV  xwJ 

uu 

u  uv  uw  — u 

C 

C  +R  C  Z 

vu 

W  V  VW  —V 

C 

C  C  +R  Z 

wu 

WV  WW  w  — w 

C  = 

C 

-  [C  C 

C 

]  C  +R  C  C  _1 

c 

XX . uvw 

XX 

XU  XV 

xw 

uu  u  uv  uw 

ux 

C  C  +R  C 

C 

VU  VV  V  VW 

vx 

C  C  C  +R 

C 

WU  WV  WW  w 

wx 

For  such  situations  the  GEOFAST  algorithm  of  Ref.  6  is 
appropriate . 
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APPENDIX  B 

ASTRONOMIC  LATITUDE  DETERMINATION 
WITH  GYROSCOPIC  SYSTEMS 
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ASTRONOMIC  LATITUDE  DETERMINATION 
WITH  GYROSCOPIC  SYSTEMS 

ADDRESSED  AS  ONE  OF  SEVERAL  ACTION  ITEMS  FROM  ETL/TASC  MEETING  OF  19,20  AUGUST  1982. 


GYRO  SCALE  FACTOR  ERROR 


FIELD  CALIBRATION  TECHNIQUES  WOULD  NEED  TO  BE  DEVELOPED; 
SEVERAL  APPROACHES  OFFER  PROMISE 


GYRO  DRIFT 


GYRO  SCALE  FACTOR  ERROR 


NET  ASSESSMENT  OF  GYRO  TECHNIQUE 
TO  DETERMINE  LATITUDE 
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APPENDIX  C 
TRANSLOCATION 
OF  DEFLECTION  MEASUREMENTS 


EXTRAPOLATION  FROM  MASTER  STATION 
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FOR  APPLICATIONS  INVOLVING  DOV  TRANSLOCATION 
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